The impact of cosmic neutrinos on the gravitational-wave background 
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We obtain the equation governing the evolution of the cosmological gravitational-wave back- 
ground, accounting for the presence of cosmic neutrinos, up to second order in perturbation theory. 
In particular, we focus on the epoch during radiation dominance, after neutrino decoupling, when 
neutrinos yield a relevant contribution to the total energy density and behave as coUisionless ultra- 
relativistic particles. Besides recovering the standard damping effect due to neutrinos, a new source 
term for gravitational waves is shown to arise from the neutrino anisotropic stress tensor. The 
importance of such a source term, so far completely disregarded in the literature, is related to the 
high velocity dispersion of neutrinos in the considered epoch; its computation requires solving the 
full second-order Boltzmann equation for coUisionless neutrinos. 
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I. INTRODUCTION 

An important discriminator among different models for 
the generation of the primordial density perturbations is 
the level of the gravitational-wave background predicted 
by these models. For example, within the inflationary 
scenario the tensor (gravitational-wave) amplitude gen- 
erated by tiny initial quantum fluctuations during the ac- 
celerated inflationary expansion of the universe depends 
on the energy scale at which this inflationary period took 
place, and it can widely vary among different inflation- 
ary models P, . On the other hand, some alternative 
scenarios, such as the curvaton model, typically predict 
an amplitude of primordial tensor modes that is far too 
small to be ever detectable by future satellite experiments 
aimed at observing the B-mode of the Cosmic Microwave 
Background polarization. 

There is however another background of stochastic 
gravitational waves of cosmological origin. Gravitational 
waves (as well as vector modes) are inevitably generated 
at second order in perturbation theory by scalar density 
perturbations [1, 0, H, [E ■ This is due to the fact that 
the non-linear evolution always involves quadratic source 
terms for tensor (and vector) perturbation modes made 
of linear scalar (density) perturbations. 



Since the level of density perturbations is well de- 
termined by CMB anisotropy measurements and large- 
scale structure observation [1, Q, we know that these 
secondary vector and tensor modes (produced after the 
primordial curvature perturbations have been generated) 
must exist and their amplitude must have a one-to-one 
relation with the level of density perturbations. In this 
sense, the scalar-induced contribution can be computed 
directly from the observed density perturbations and gen- 
eral relativity, and is independent of the specific cosmo- 
logical model for generating the perturbations.^ 

Such a background of gravitational waves could be in- 
teresting in relation to future high-sensitivity CMB po- 
larization experiments or for small-scale direct detectors, 
such as the space-based laser interferometer Big Bang 
Observer (BBO) and the Deci-hertz Interferometer Grav- 
itational Wave Observator y (' DECIGO) operating in the 
frequency range 0.1 - IHz [111 ] with an improved sensitiv- 
ity (in terms of the closure energy density of gravitational 
waves, ^GW ^ 10"^'' — 10~^^). In particular, in Ref. [12 
the effects of secondary tensor and vector modes on the 
large scale CMB polarization have been computed, show- 
ing that they dominate over the primordial gravity-wave 
background if the tensor-to-scalar perturbation ratio on 
large scales is r < 10^^. More recently, Ref. [l^ com- 
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^ See however Ref. [loll , where in the context of the curvaton mech- 
anism, second-order gravitational waves can be produced when 
the perturbations are still of isocurvature nature, thus resulting 
to be strongly model dependent. 
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puted the power-spectrum of the secondary tensors ac- 
counting for their evolution during the radiation domi- 
nated epoch, to see their effects on the scales relevant for 
small-scale direct detectors, and Ref. [TJ] extended this 
analysis by accounting for a more detailed study of the 
transfer function for the secondary tensor modes. 

In this paper we consistently account for the presence 
of cosmic neutrinos to analyze their impact on the evolu- 
tion of the second-order gravitational-wave background. 
At linear order it has been shown that there is a damping 
effect due to the anisotropic stress of free-streaming neu- 
trinos that strongly affects the primordial gravitational- 
wave background on those wavelengths which enter the 
horizon during the radiation dominated epoch (at the 
level of 30%) ji5|,lil,[i3,[li,[i9|,|20j,l2l| (see also Ref. [13]). 
At second order, along with the analogous damping ef- 
fect, we find that free-streaming neutrinos are an impor- 
tant source for the second-order gravitational-wave back- 
ground during the radiation-dominated epoch. We find 
completely new source terms, arising because of the fact 
that neutrinos give a relevant contribution to the total en- 
ergy density during this epoch and they behave as ultra- 
relativistic collisionless particles after their decoupling: 
their high velocity dispersion acts as an extra source for 
the second-order gravitational waves. To compute such a 
contribution we evaluate the second-order tensor part of 
the neutrinos' anisotropic stress tensor, that has been ne- 
glected so far. This is achieved by computing and solving 
the Boltzmann equation for neutrinos. Approximating 
the neutrino contribution as a perfect fluid of relativis- 
tic particles during the radiation era leads to a serious 
underestimate of their role. Let us stress that the new 
contribution is at least of the same order of magnitude 
as that computed by adopting a fluid treatment in the 
source of the scalar-induced gravitational waves. More- 
over it has a clear physical interpretation. It arises in the 
Boltzmann equation from a "lensing" effect of the neu- 
trinos as they travel through the inhomogeneities of the 
gravitational potential. 

The paper is organized as follows. In Sees. Ill Al and 
IIIBI we derive the evolution equation for the tensor 
(gravitational-wave modes) at second order, accounting 
for photons and neutrinos. In Sec. IIII Al we present 
the Boltzmann equation for neutrinos approximated as 
being collisionless massless particles (see Appendix A, 
B and C for details about the Boltzmann equation for 
massive neutrinos) and give an integral solution for it. 
Sec. IIII C] contains the computation needed to determine 
the tensor part of the second-order anisotropic stress ten- 
sor of neutrinos, which leads to one of our main results, 
Eq. ([5^ . Finally, in Sec.|lV]we derive the photon contri- 
bution, consistently accounting for the presence of neu- 
trinos, which leads also to a new source of gravitational 
waves. In Sec. |V] we present the summary and our main 
conclusions. 



II. SECOND-ORDER GRAVITATIONAL WAVES 

A. Metric perturbations in the Poisson gauge 

The second-order metric perturbations around a flat 
Friedmann- Robertson- Walker (FRW) background can be 
described by the line-element in the Poisson gauge 

ds^ = a2(T) [-e^^dr^ -|- 2uj,dx'dT + (e"^*^^^. 

-I- Xi])dx'dx^] . (1) 

In this gauge one scalar degree of freedom is eliminated 
from goi and one scalar and two vector degrees of free- 
dom are removed from gij. As usual a(r) is the scale 
factor and r is the conformal time. The functions and 
^' are scalar functions which correspond to the Newto- 
nian potential and to the spatial curvature perturbations, 
respectively. Within second-order perturbation theory, 
they consist in the sum of a linear and a second-order 
term, such that $ and "if can be written as 

$ = $(1) +$(2)^2 and = ^'(1) +<If(2)/2. (2) 

Since the choice of the exponentials greatly helps in sim- 
plifying the computation of many expressions, they will 
be kept where it is convenient. It is worth remarking here 
that all the equations in the following where the expo- 
nential show up are meant to be second order equations, 
therefore the exponentials are to be thought as implicitly 
truncated up to second order in all these expressions e.g. 
e2* ~ l + 2$(i) + $(2) + 2($(i))2. 

The remaining functions that appear in Eq. ^ account 
for second-order vector (uji) and tensor (xij) modes. Ten- 
sor perturbations are traceless and transverse: xl — Oj 
diX^^ — and vectors have vanishing spatial divergence: 
d^LUi = 0. Linear vector modes have been neglected as 
they are not produced by standard mechanisms, such as 
inflation, that generate cosmological perturbations [23j . 
dl]. As discussed, for example in Refs. [1^ and [13], in- 
deed linear vector modes have decreasing amplitudes and 
they are not generated in the presence of scalar fields, 
while the first-order tensor part gives a negligible con- 
tribution to second-order perturbations. Second-order 
vector and tensor modes however must be taken into ac- 
count, even if they were initially zero. This is because 
scalar, vector and tensor modes are dynamically cou- 
pled at this stage and second-order vectors and tensors 
are generated by first-order scalar mode-mode coupling. 
First-order perturbations behave as a source for the in- 
trinsically second-order fluctuations ^ . 

Since our main task is to provide the second-order 
Einstein's equations that describe the evolution of ten- 
sor modes, we are interested in the spatial components 
of both the Einstein and the energy-momentum tensors. 
Here we find that, using the Christoffel symbols obtained 
in Appendix and accounting only for the terms up to 
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second order, the spatial Einstein tensor reads 



-2* / n_,2 



n-2 — 2^-'$' - 3(^'y + (3) 

a 



2H($' + 2*') + 2^" I + e^* f dk^d''^ 



„2* 



,2 y^x)+^x) , 



where Ti. = a' /a, and a prime denotes differentiation 
w.r.t. conformal time. 



B. Second-order gravitational- wave evolution 
equation during the radiation-dominated era 

In Fourier space the equation which describes the evo- 
lution of second-order gravitational waves (GW) can be 
put in the form: 



Xk,x + 2WXfe,A + k'^Xk,\ = IGTrGa^S, 



(4) 



where the subscript A refers to the two possible polariza- 
tion states of a gravitational wave. Each mode Xk is in 
fact transverse with respect to the direction along which 
it propagates and, for a mode traveling in the z direction, 
Xij can be written as: 



Xij 




(5) 



The two degrees of freedom account for the two polariza- 
tion states x+ and Xx. 

The source term Sk,\ for GW in the radiation era con- 
sists in the sum of three different parts: Ek.\ that comes 
from the Einstein tensor, that comes from the neu- 

(7) 

trino anisotropic stress tensor term and Tlj, ^ that ac- 
counts for the photon contribution. We then have: 



Sk,x = Ek,x + n 



fe,A 



n 



(7) 

fe,A 



(6) 



erator V^y , so that 



2. 



(7) 



Here T''^ contains both the neutrino and photon con- 
tributions 

The definition of such an operator is given in Ref. [6| 



T) is _ 'pi'ps 'pi'ps 



(9) 



where 

vi^5)-{^^Y^ d'dj. (10) 

Moving to Fourier space the two-indices operator T'* 
reads 

V]^[b)~k%), (11) 
and then the operator to apply is 



(12) 



The contribution to the source term coming from the 
Einstein tensor can then be written as: 



,2 



(13) 



To get Ek. X we now have to move to Fourier space and 
project along one of the two polarization states A = -I- , x . 

If we take an orthonormal basis made up by the three 
unit vectors e, e and k, the two polarization tensors are 
defined as follows: 

4(k) = -L[e,(k)e,(k) - e-.(k)e-,(k)] , (14) 

£,^.(k) = -L[e,(k)e-,(k) + e-.(k)e,(k)] . 

Remembering that a gravitational wave is transverse 
with respect to the direction k along which it propagates 
we find that : 

(fkld^k2 c., , , , ^ 

-3— d(ki ka - k) 



Ei = e'^E' ---2 



^ " ' '2Tr) 

A:ifc2$fei$fe2 - 2kl^k^'^k2 + 2fc^*fci^'fc2 

/Cife^fci^'fca - ^lfc2*fei^'fc2 + fclfe^fci^fesl • (15) 



X e{ kijk^ 



In making the source term Sk^x explicit, the first step 
is to extract the tensor part of the Einstein and energy- 
momentum tensors to get the corresponding transverse 
traceless component. 

This can be done by making use of the projection op- 



^ In the following, where not necessary, we will omit the TT sub- 
script. 

^ The indices of the 3-D unit vectors are lowered and raised by 
using 5ij so that S^e^Cj = 1 = (5*e'ej. 
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The explicit form for the product sl^kijkl^ is given by 
the analogous of eq. 

The next sections are devoted to the computation of 
the neutrino source term since the term due to 

the photons is already known and we will recall it later. 
Notice however that Sec. IIVI presents a non trivial case 
accounting for both neutrinos and photons: here in fact 
the contribution of the first-order coUisionless neutrinos 
is fully accounted for the first time when explicitly cal- 
culating the second-order photon quadrupole in the tight 
coupled limit (Eq. [59]). 

We will proceed by steps. First we solve the second- 
order Boltzmann equation for the neutrinos in order to 
give the expression for the neutrino energy momentum 
tensor. Then we will give the expression for the spatial 
components of the energy momentum tensor of neutrinos 
in terms of the their distribution function perturbed up 
to second order, and finally we can extract the transverse 
and tracelees part of it. 



III. NEUTRINOS AND SECOND-ORDER 
TENSOR MODES 

A. Solution of the second-order Boltzmann 
equation for neutrinos 

The Boltzmann equation up to second-order for decou- 
pled neutrinos can be written as 



dF^ dF^ dx' dF^ dE dF^ dn^ 
dr dx^ dr dE dr d-n} dr 



0, (16) 



where r is the conformal time, E — ^ rn^ -V , = 
gijP^P^ is the squared neutrinos 3-momentum and is 
the neutrino 4- momentum defined as — dx^ /d\. Here 
A parametrizes the particle's path and x^ — {t, x) repre- 
sents a space-time point. The unit vector n, with com- 
ponents n*, represents the neutrino momentum direction 
and it is therefore such that: p = pn and rfn^Sij = 1. 
Equation pB]) refers to the general case of massive neu- 
trinos igf,aP''P" = -ml). 

The neutrino distribution function F^, at this stage will 
consist in an unperturbed, a first-order and a second- 
order term, such that we can put it in the form 



(2) 



F, = F,{E, t) + F(,^\E, T,x\n') + ^{E, T,x\n'). 

(17) 

The explicit forms for the different contributions in 
Eq. (flB)) are listed in Appendix [Cl 

Since we are interested in evaluating the neutrino con- 
tribution during the radiation-dominated epoch, we can 
specialize the second-order Boltzmann equation for ul- 
trarelativistic (massless) neutrinos. Replacing E with 
the comoving 3-momentum q as one of the independent 
variables in Eq. (|16p such that p = E and q = ap, the 



relevant equation takes the form: 



IdF^ 



(2) 



(2) 



2 di 



2 dx' 



-n 



dFk 



(1) 



dx^ 



.„*($(i) +$(1)) 



dF 



vi/d)') 



„J(l$(2)+ci,(l)(vI/(l)+$(l))) 
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Going to Fourier space it can be put in the form: 



1 



P, 



(2)' 



zfc/iP^'^ -Gk(T)-Tk(T), 

/Lt = k • n, 



where 

and 
Gk(r) 



dF 



(2) 



(2)' 



(18) 

(19) 
(20) 

(21) 



d^kid''k2 ^on. , L , ,t,(i)^ 



(27r) 



- iVi ■ n- 



9P 



(1) 



dri 



9P, 



(1) 



9- 



^(-lfc2/i$i'^+* 



,(1). 

- ki > 



Here Vn accounts for the second-order vector pertur- 
bations modes, see Eq. (|C17p . We will keep it in this 
implicit form, since we already know that second-order 
vectors do not take part to the tensor contribution we 
are interested in. 

The term Pk represents the "pure" tensor contribution 
and is defined by 



dF 



(22) 



A=x,+ 

The Fourier expansion for the tensor modes reads 

d^k 

A=x,+ 



(2^; 



3Xk,A(r)e^^=^4 (23) 



where e^- are the polarization tensors defined in Eq. p4|) 
with A = X , + accounting for the two possible polariza- 
tion states of a gravitational wave. 

At this point, following the standard procedure [l^, 
[28| , we can write down an integral solution of the second- 
order Boltzmann equation as 



rf^'e»fcMr'-r)[G,(/)_Tfc(r')]. (24) 
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As usual, Tdec refers to the neutrino decoupling conformal 
time 

For F^'^ we take the formal solution of the first order 
neutrino Boltzmann equation integrated by parts: 



(25) 



where /^i = ki • n. In order not to weigh further the 
notation, will be explicitly inserted only later, when 
strictly necessary. 



B. The second-order neutrino energy- momentum 
tensor 



The contribution to the energy momentum tensor of a 
given species "z" is 



g J (27r)3 PO 



(26) 



where P'^ = dx^ / d\ 4-momentum of the particle and Pi 
is the distribution fucntion of the given species. 

If we now want an expression for the second-order spa- 
tial component of the neutrino energy momentum tensor, 
we find that it will consist of the sum of four parts. There 
are in fact four terms which contain the perturbations: 
the product of the 3- momentum P'P^, the component 
, the distribution function F^, and the determinant of 
the metric ([I]), 5. As far as the first and the second term 
are concerned we must remember that 



= ^e* ( 1 - ^Xmin^n 



1 



(27) 



while for the distribution function we use the decompo- 
sition in Eq. (jl7p . By performing the variable change 
Pj Qj in order to make all the perturbations explicit 
in eq. (|26p and by combining all the terms, we find that 
the neutrino energy momentum tensor at second order in 
perturbation theory reads: 



-4 

a 



dh 



(2^) 



(2) 



(31) 



Similar to the linear case, vectors and tensors do not show 
up because of the angular integration. 

With this results we are able to express the spatial 
component of the second-order energy-momentum tensor 
Eq. (|3ip in Fourier space: 



(2) 



2a' 



d\ 



-qn n, 



(2^) 

rfTV'=^(^'-^)[Gfe(T')-Pfe(T')] 



(32) 



C. The neutrino contribution to the source term. 

In order to get the transverse traceless part of the neu- 
trino energy-momentum tensor we now have to make use 
of the operator defined in Eq. (|12p. Before proceeding 
notice that we are interested only in the tensor contri- 
bution to the energy momentum tensor. Therefore it is 
very useful to keep in mind the decomposition of the dis- 
tibution function into its scalar, vector and tensor parts, 
according to the splitting of Ref . [22] . For example the 
tensor part is given by 



5F 



A 



/A(k,r,g,n)e,^ 



(33) 



2uj, 



- UJiTl 



(28) 



(29) 



where we define = a'giji 
momentum q = gn of magnitude q and direction n*, see 
the notations of Ref. [25] and [26]. The overline refers 
to unperturbed quantities and we have q^ = a^^ qn'' , 
Tjj = aqrij and if = q/a. 

The determinant of the metric, g, up to second order 
is such that: 



a-4e3*-* 



(30) 



This greatly helps in simplifying all the expressions: it 
is telling us that the only tensor contributions to the 
energy momentum tensor are those that can be built out 
of the product of two momentum direction n*. Looking 
back to Eq. (pij) and Eq. we see that, apart from 
the straightforward term depending on the gravitational 
waves Xijj there is just one term of this type, namely 



ap^ 



(1) 



(34) 



Therefore, from now on we will focus just on this term in 
Eq. and we will drop all the others since they cor- 

respond to scalar and vector modes. ^ It is interesting to 
notice that the term (I34p has a clear and simple physical 



* Notice that we are assuming that at a time right before T^g,^ the 
distribution of the neutrinos is the homogeneous one since they 
are still in thermal equilibrium. 



^ In fact we have explicitly verified that all the remaining terms 
vanish once the energy momentum tensor is projected along 
the two polarization tensors l|14|l and the angular integration 
in Eq. II32I I is performed. 
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interpretation. It arises in the Boltzmann equation from 
a "lensing" effect of the neutrinos as they travel through 
the inhomogeneities of the gravitational potential. In the 
Boltzmann equation (|16p it derives from Eq. (jClip . which 
describes how the neutrinos momentum direction changes 
in time due to the potential wells they pass trough. 

Let us now continue our computation and apply the 
projection operator (fT^ to Eq. ([5^ . keeping only the 
term (p4)) in Gk('r') and the "pure" tensor contribution 
Tw{t'). We see that this operator will act on the product 

of the two direction unit vectors n^Ug contained in Ts^'^\ 
Since S^^.n^rig — 1 and the product k-p defines the cosine 
/i between the direction along which the perturbation 
propagates and the neutrino momentum, this means that 
(niOk"^ will contain the term: 



n'nj-n{kjn' + k'nj) + -k'kj{l + ^i^)^-d){l-fi^). (35) 



As already done in Sec. IIIB| we can now choose one 
of the two polarization states A = + , x and project the 
transverse traceless part we are interested in into the cor- 
responding polarization tensors . Since a gravitational 

wave is transverse with respect to the direction k along 
which it propagates, we then have 



1 



Si^n'rij = ^(l-/)((5+^cos2v3n + (5''-^sin2v?„), (36) 
V 2 

where the angle ip^ is the azimuthal angle of the neutrino 
momentum direction n in the orthonormal basis e, e and 
k. In Fourier space the term p4|) evaluated at r' becomes 
(for simplicity we omit the convolution integral) 

- (ki.n)(k2-n)($k.(r') + *k.(r'))g|^ (37) 
X r dT"e**'-i'^i(^"-^')(r"-r')«,(T") + <,(T")), 



where we have written explicitly the term 9-Fki / dn' as 



dn'- 



dr'zA:i(r-r')K(r') + *',,(r' 



X e 



ifci/^i(r'-r) 



(38) 

using the solution of Eq. (PS)) . It proves convenient to 
take the tensor part of (ki • n)(k2 • n) which reads 

(fc2ifci™e+"")e+ n'^n^ + (fc2ifci™£^'™)e.^.n'-n^ . (39) 

Notice that this step js equivalent to follow the decompo- 
sition ([33]) of Rcf. 22] which allow to isolate the tensor 
contributions to the distribution function.^ 



With these results we are now able to write the expres- 



sion for = ei^"'II'^j^'' at second order in the perturba- 



tions 



n 



k.A 



= -2a" 



d^ki Sk2 



(27r)3 



<5(ki + k2 - k) 



d^g 2 dF 



X [{k2ikime+'n4sn''n' + (A:2iA;i™e^'")e,^,n'-n^] 
X A^2{t') r dr"e*'=i'^i(^"-^')(r"-r')A;,,(r") 



+ a g. 



(27r)3 4 dq 



dr' 



(40) 



where A^ir) = $k(T) + ^'k(T). 

At this point we want to solve the angular and the mo- 
mentum dependence of the transverse, traceless neutrino 
energy momentum tensor Tl'^'^ ■ This implies solving the 
d^q integral in the previous expression. 

First of all notice that the pure tensor part in the last 
line of Eq. (|^ represents the second-order analogue of 
the damping effect found and discussed in Refs. [19, 
[13, [H [H, 0, [m, ill . We are then able to deal with it 
(see, for example, Refs. [IMIllI, 113) ^tnd we find 



dqq^ dF r , , ^, 1 , . , , 10 . , , 3 . , 

[T^r dq Ju^^^ lb 7 7 

(41) 



which becomes 
-8p.(r) 



(42) 

where we have used Eqs. (|50p - ([5T|) . the variable change 
T u — fcr, s — u — U, with U — /ct', and the derivative 
is with respect to U . This means that also at second 
order gravitational waves are damped by neutrino free 
streaming, as expected on general grounds. Here py is the 
unperturbed neutrino energy density given in Eq. (|5D1) . 

We now have to deal with the reminder of Eq. (|40p . 
giving an additional source term that represents a com- 
pletely new result. 



1. Angular integration 

We notice that there is a /i-dependence hidden in the 
exponential 



Also in this case we have verified that the scalar and vector com- 
ponents of (ki • n)(k2 ■ n) give a vanishing contribution once the 
energy momentum tensor is projected along the two polarization 
tensors II14I I and the angular integration in Eq. I I32I I is performed. 



^z'(2/-M)j,(s')^/(Mi) (43) 

+ 00 +1 

4^E*'J'(^') E «;™e''"'^"Pi™(A^)F,:„(ki), 



1^0 



7n— — l 
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where ai, 



(2i + l)(i-m)! 
(i+m)! 



and s' = fci(T' - r) = {Ui - 



Ml). For the integration over the angle ip^ we need to take 
into account the product (ef"' n^nj){e^gn^n^) in Eq. ([40]) 
and, using Eq. (|36p . we find 

d(/9„e""'^"(£,^-'' n'n,)(£^;n'n^) - |(1 - ^2)2^aa'^^^^^ ^ 

(44) 

which means that the cross terms vanish, while the 
squared ones select m — 0. When m = the P;,„ corre- 
spond to the Legendre polynomials Pi and we have: 

-t-oo 

v^j2ji{s'yv2r+TYi*{ki)Pi{^i). (45) 

1=0 

Therefore by expanding e*'^''*-'^ as in Eq. ([43]) we are 
led to an integration over /i of the following quantity 



J -I 111 



(46) 



where 



Ai - zV4^(2? + l)r,5(ki)j,[fci(r"-r')], 
B,, = i''(2Z' + l)>[fc(^'-^)]- (47) 

Such an integral is easily performed using the formulae 
of Appendix [Pl and Eq. (|46|) becomes 



E 



1 



I 



21 + 1 



(48) 



where 



A 



(1) 



Z(/-l) 



(2;-3)(2Z- 1) 



Al-2 + 



(2; + 1)(2Z- 1) 



A, 



{2l + l){2l + i) 
{l + 2){l + l) 



(2Z + 3)(2Z + 5) 



A 



1+2 ■ 



(49) 



and similar for 



(1) 
; ■ 



3. Final expression for the neutrino contribution 

Collecting the previous results we arrive at the contri- 
bution of streaming neutrinos to the evolution of gravi- 
tational waves: 



E 



(2^)^ 

dr'(<i>k,(T') + *k.(T')) 



(1) 



(1) 



X / dr"(r"-r')«,(r") + <,(r")) 

/ 1 10 3 



(52) 



where s — k{T — t'), and the functions Ai^ Bi and 

are defined in Eq. (gT]) and Eq. 

This equation is clearly telling us that coUisionless free 
streaming neutrinos contribute with new terms to the 
source of second-order gravitational waves, with respect 
to the fluid treatment adopted in the literature so far, 
where the tensor part of n'^''^ at second order has never 
been taken into account. 

Even with a qualitative approach, we can state that 
if low / contributions (up to ? = 2) can have a corre- 
spondence with respect to a source where neutrinos and 
photons are treated as a single-fluid radiation, higher 
multipoles surely do not. These in fact come from the 
high neutrino velocity dispersion and can be found only 
if neutrinos are treated as coUisionless particles. 

Moreover, during the radiation-dominated epoch, the 
neutrino fraction fu(j) = Q.^/Vlif is not negligible. Since 

all the terms in li!"^^ are multiplied by p^, they are 
therefore non-negligible at that time. 



IV. THE PHOTON TENSOR QUADRUPOLE 



2. Integration over the comoving momentum q. 

We can treat the momentum integration over q in 
Eq. (j40p independently from the angular part. There 
is just one type of q-dependence, namely q^dF /dq. Re- 
member that F is the unperturbed neutrino distribution 
function given by the Fermi-Dirac distribution and 



-4 3* 
= " 2^ 



Fq^dq^l—Tt, 
^ ^ 8 30 " 



(50) 



is the unperturbed neutrino energy density. Hence, inte- 
grating by parts when necessary, we find (the 7r/2 factor 
comes from Eq. dM])) 



dq dF 
(2^ 9^ 



1_ 



(51) 



In order to complete the expression of the source 
term for gravitational waves during the radiation era in 
Eq. ([¥]), we now have to add the photon contribution. 
Following [2^, the second-order photon quadrupole in 
the tight coupling limit is given by: 



n 



(2)U 



(53) 



To get the photon contribution to the second-order 
gravitational waves source we then have to extract the 
transverse and traceless component. 

Since the operator to apply is such that V^fSl — 0, we 
will find that the tensor quadrupole takes the form: 



«?) 



TT 



-V 1,(1)'- 1,(1) 
3 i-J 7 7 s 



(54) 
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By making use of the first-order space-time component 
of the Einstein equations (see Appendix 



1 



WttG , 



we can define the first-order photon velocity as: 
3 



WttGo? p. 



(55) 



(56) 



where v^'^^ represents the first-order neutrino velocity. 
During the radiation-dominated epoch neutrinos are still 
relativistic so that, in Fourier space, the velocity can be 
written as: 



1 



/ (0^"^'^^' = -l2niM^k\ (57) 



Pu 



p. 



The term Mi refers to the neutrino dipole and is the 
neutrino pressure. By making use of the formal solution 
P5|) of the first-order neutrino Boltzmann equation, the 
dipole can be expressed explicitly in terms of the pertur- 
bations to give: 



idlnF, 



2 d\nq 
**fciJi(s) 



dr' 



^$',^(jo(s)-2j2(s)) 



(58) 



the evolution of the gravitational wave background; it 
provides, first of all, a complete study of the Boltzmann 
equations for neutrinos at second order and the expres- 
sion for the second-order anisotropic stress tensor. 

Free-streaming neutrinos are an important source of 
second-order gravitational waves during the radiation- 
dominated epoch. Along with the fact that neutrinos 
yield a relevant contribution to the total energy during 
this epoch, this is due to the large neutrino velocity dis- 
persion and it emerges from the calculations performed in 
this paper when assuming that neutrinos are coUisionlcss 
particles, as it is the case after their decoupling. 

The fluid treatment adopted so far for describing 
neutrinos turns out therefore to be a poor approxima- 
tion that leads, in particular, to underestimating the 
role of neutrino free-streaming as a source of gravita- 
tional waves. In particular, in this paper we have made 
the first full consistent computation of the second-order 
tensor part of the neutrino anisotropic stress tensor, 
Eq. ([5^ . This has been achieved by computing and 
solving the second-order Boltzmann equation for the 
neutrino distribution function. Besides recovering the 
second-order counterpart of the damping effect studied 
in Ref. [H, M, [13, M, M, S mi, m, Eq. ^ repre- 
sents a completely new source term for the evolution of 
gravitational waves. 



Projecting along one of the polarization state (-I-), we 
have that, in Fourier space, the photon contribution we 
are searching for takes the form: 



n 



(7)A _ 



(27r)3 
3fcifc 



(5(ki+k2-k)£^'^fcljfc^ 



32(7rGa2p ) 



+ 32-12 (^^j M(fci)M(fc2) 
Pi 



k2 



(59) 



QPu 



^fciAAi(fc2)(*'fc^ 



Ga-^p 



It is worth stressing that this expression contains at 
most multipoles with I = 2. Therefore, if we had 
used a fluid treatment for the neutrinos as well (e.g. 

cx vi^^\i^^^), we would have found in the second- 
order neutrino contributions terms with I not higher then 
1 = 2. This is clearly in contrast with the neutrino source 
term Eq. ((52)) . 



V. CONCLUSIONS 

This paper represents the flrst step towards the quan- 
titative evaluation of the impact of cosmic neutrinos on 
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APPENDIX A: THE SECOND-ORDER EINSTEIN 
TENSOR 

In this appendix we provides the deflnitions for the 
connection coefficients and the expression of the second- 
order Einstein tensor for the metric ([1]): 



ds^ = a^(r) [-e^^dr^ + 2ujidx'dT + {e-'^"^ Si^ + Xzj)dx'dx^ 



-2* i 



The space-time metric g^^^, has signature (— ,- 
The connection coefficients are defined as 



^01 — off 



2 ^ \dxP dx-r dxP 



(Al) 



Greek indices {a, P, p,v, ....) run from to 3, while 
latin indices {a,b, . . . j,k, . . . , m, n, . . . ) run from 1 to 
3. In particular, their explicit expression with our metric 
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reads: 



(A2) 



00 



ax,; 



1 / dOJn 



+ 



2 V dx^ dxi 

2 

1 



+ e 



-2*-2* 



- jk 



a* 9* 
-Huj'Sik - —5\ - —6\ + —5jk 



2\dx'' dx:> 



dx'^ ■' dx^ 



dxi 



The Einstein equations are written as G^^ = 

SttG^T^^, where Gn is the usual Newtonian gravitational 
constant, G^ji/ = R)j,i, — \g^.uR is the Einstein tensor and 
Tpjy is the Energy-momentum tensor. 

The Ricci tensor R^i, is a contraction of the Riemann 
tensor, R^i, = R^^^iu '^'^'^ ™^ terms of the connection co- 
efficient it is given by 

p p) pa Pi -pet I pa p<7 pa ptr ( /^"i^ 

J^ixv — '-'a'- '-'IJ-'- va ^ '- era'- tiv av na ■ \-^'^) 

The Ricci scalar is the trace of the Ricci tensor, R = i?^., . 
The components of Einstein's tensor up to second-order 
read 

G'o = - 



„2<I' 



(S^f^^* - 2V2^')] , 



2— [d'^' + {H- 5*$] - 



+ (m^- 2- 



(A4) 



(A5) 



G\ = 



-2$ 



-2—- 2*'$' - 3(*')^ + (A6) 



+ 2H{^' + 2*') + 2*"^ -I- e^* {dk^d^^ + 



o2* 



The exponentials are maintained since they help in sim- 
plifying lots of calculation, however notice that in all 
these expressions they are implicitly truncated up to sec- 
ond order. 



APPENDIX B: FIRST-ORDER 
PERTURBATIONS OF EINSTEIN EQUATIONS 
FOR PHOTONS AND NEUTRINOS 

Carrying out the calculation, the first-order Einstein's 
equations, expressed in Fourier space in terms of the per- 
turbations $ and ^, take the form: 



-fc2*-3^ f * - $^ ) = 16nGa^[p^e^o + P,Q.o] (Bl) 



A;2($ - ^r) = -327rGa2p^e^2. 



(B2) 



Here Q^q and O^^q ^-rc, respectively, the photon and 
the neutrino monopole contribution, while 0j,2 refers to 
the neutrino scalar quadrupole. 

Remember that the first-order perturbation to the neu- 
trino distribution function is defined by 
where 



= — Si ^'^ 

omp 



(B3) 



and Qi,{t, x, n) = Similarly it happens to the photon 
contribution. 

In general, the Zth multipole of the temperature field 
Q can be defined as: 



{-^y j-1 2 

where Pi (/x) are the Legendre polynomials. 
The neutrino monopole is then defined by: 

while the quadrupole, corresponding to i = 2, is 



(B4) 



e.2 = -/^f P2(M)e.. 



(B5) 



(B6) 



APPENDIX C: THE SECOND-ORDER 
NEUTRINO BOLTZMANN EQUATION 

In this appendix we provide the explicit form for the 
different contributions in the neutrino Boltzmann equa- 

+lr.Ti -1- ^^i' rf^i J- ^^k dE I dF^ dn' _ n . 

I'lOU -t- -g^ ' dr ^ dE ' dT ^ 'dn^ ' dr ~ " • 



dx^ 



dx^ d\ — P' 



dr dX dr — P°- 

For this term we have: 



dr E 



i P 
■ OJi w — 

^ E 



(CI) 
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dE 
dr ■ 

Deriving an expression for this term is more lengthy 
since it involves the use of the geodesic equation. 

We can start noting that, using Eq. (|29p. we obtain: 



dr 



E 



d<P dB /I dE 1 da 

hS 

dr dr \ E dr a dr 



, (C2) 



Carrying on the calculation in the same way as we 
did for dE/dr, at the end we recover the expression 



dP' _ p 
dr a 



dr dr p \ E 

E'f')] . (CIO) 



where B 



1 + %UJ^n' 



To simplify a bit the notation, from now on we will 
set 



, _ dg _ dg 

g=- and 9, = g-, V.g 



(C3) 



If we now make the total derivatives explicit, ac- 
cording to the fact that 



d7 = ^+*-P^^ 



and 
dB_ 

we find 

1 dP^ 
Elh' 



p idui 



E 



J po J 



(C4) 



(C5) 



dE 
d7 



1 da 
a dr 



P^ , 



■^1,3 



By making use of the geodesic equation, we can 
express the time component as a sum of three terms 



The term we are looking for has the form 



dn^ 



:^$.'_:^-jf^» + „^nJ"-$ , + k- (CU) 

„ 'J p ' ' ' 



P 



P 



P 



It is worth noticing that, since dF^/dn^ is already 
a first-order term, we must consider this equation 
just up to first-order. The scalars that appear here 
are therefore the linear components of the pertur- 
bations. 

The purpose of this Appendix is to find the second-order 
Boltzmann equation for (decoupled) neutrinos. It is then 
useful to use the relations in Eqs. (fT7)) and ^ for making 
the second-order terms explicit. 

The second-order contributions of each part of equa- 
tion are listed below in the following expressions: 



dF^ dx^ 



) . (C6) dx^ dr 



2„.orf dx^ E '2 dx' E 

(C12) 

The energy dependence term is a bit more complicated 
since it consists of three terms such that: 



^(2) 



pa pP 

po 



pO pO 
00^ 



pa pP 
^0 ^ ^ 

V po 



(C7) 



If we now insert the connection coefficients and 
the perturbed components of the 4- momentum (I^T)) 
and (|29p . we come at last to the required term 

2 



P_ 

E 



-( 



P.YL' 

e) 2 



\E 



t„3 



We can note that the first term in the RHS of the 
latter equation is the zeroth-order time component 
of the geodesic equation. 



dF^ dE 
~dE~d^ 



dF^ /dE 
'dE [jh 

IdFi)^^ fdE 
2 dE V dr 



dF^^^ fdE 



dE \ dr 



(C13) 



(The overline refers to the zeroth-order neutrino distri- 
bution function). 
We then have 



(C8) 1 ( ^ 

E\dT 



2,^d 



p\2 1 



^(1)^(1) ^(1)^(1) 



pN^2 1 

r 



(C14) 



dn^ 



This term requires some lengthy algebra as well. 
To obtain it the spatial component of the geodetic 
equation must be used: 



dP' 
Ik 



pa p0 
a/3 po " 



(C9) 



while 



1 dE 
E'd7 



and 



E 



(1) 



E 



(C15) 
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1 dE 



Qth 



n. 



(C16) 



The term Vu accounts for the second-order vector con- 
tribution. This aheady appears in Eq. (jC8[) and is defined 
as 



Vu - ^ 



(C17) 



Finally, we use Eq. (|C11|) to deal with the dependence 
on the momentum direction. Matching together all these 
terms in Eq. (|16p. we have now all the tools needed to 
obtain the second-order Boltzmann equation. 



and similar for Bi^\ The result in Eq. ([48)1 is obtained 
using the orthogonality of the Legendre Polynomials and 
applying to Eq. the formula (here ai is a generic 
function of I) 



(D2) 



where 



APPENDIX D: SOME FORMULAE USED FOR 
THE ANGULAR INTEGRATION 



The coefficients A'j^^ and B^^'' appearing in Eq. 



are simply given in terms of Eq. ((47l) as 



l{l^l) 



(2;-3)(2Z- 1) 



Al-2 + 



{2l + l){2l + 3) 



{2l + l)(2l- 1) 



, (/ + 2)G + 1) 



(21 + 3) (21 + 5)' 



(Dl) 



I 



1(1-1) 



{2l-2,){2l- 1) 



ai-2ii-2{x) + 



{1 + ^? 



(21 + l)(2l - 1) 



aiji (x) 



(2l + l)(2l + 3) 
(l + 2)(l + l) 



(2l + 3)(2l + 5) 



ai+2ji+2(x) , 
(D3) 



which derives from the recursion relation of the Legendre 
Polynomials (see also [30l|). 
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